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The five parameter double integral /" dy exp(— p 2 y 2 ) sin(/3r + 0) times f% tix exp( — a 2 ) cos(e/3* 4- (/>) is 
evaluated in terms of Fourier transforms of exp(— ;t 2 )erfc(a*;). Some new expressions for these transforms are 
obtained. 

Key words: Definite integrals; double integrals; error functions; Fourier transforms; radiation reaction. 

1. Introduction 

The five parameter definite double integral considered in this paper is a generalization of that encountered 
in describing the radiation reaction effects on a charged particle swept over by a single plane laser pulse. The 
electric field for such a pulse can be described by the function 

E(u) =E°exp(-// 2 //> 2 )sin(///A) 

where u = ct — &'X, k being the direction normal to the plane wave. 
The integral encountered is 

M = J c£E(£) J E(r))-E(r))dT} 

which in the notation of this paper may be expressed as 

M = V 2 D 2 (E° -E )E°[/( 1/V2, \/2/c, 0; 0, 0) - /(1/V2, V2/c 0; 2. 0)] 
and evaluated to yield 1 

M = V 2 D 2 (E° •E°)E°[exp(-K 2 )erfi(\/2K/\/3) - Vi exp(-3K 2 )erfi(2\/2K/\/3)] 

where k= D/{2\). 

The five parameter definite double integral containing Gaussian and trigometric factors, discussed in this 
paper is 



/(/>.£. 0; e,0)=| dyexp{-p 2 y 2 )sin(Py + 6) dxexp(-x 2 )cos(epx + (f>) 



(1) 



where p > 0, ft > 0, and e ^ 0. A somewhat related indefinite double integral, containing only Gaussian 
factors, has previously been discussed by Rosser [3]. It will be shown that/(/>, /3, 0; e, (/)) may be evaluated 
in terms of Fourier transforms of Gaussian weighted erfc functions. The relevant Fourier sine transforms are 



* Invited paper. 

1 The normalization adopted for the erf, erfc, and erfi functions in this paper is the same as that employed in references [1] and [2]. Figures in brackets 
indicate literature references at the end of this paper. 
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available in closed form. Closed form expressions of such Fourier cosine transforms also exist for special 
cases, and integral forms convenient for numerical evaluation may be developed for the general case. The 
necessary Fourier transforms are discussed in the next section. 

2. Fourier Transforms of exp(— x 2 )erfc(ax) 

The Fourier sine transform, for a. > 0, 

F s (2jS,a) = I sin(2j3*)exp(-^)erfc(a*) dx (2) 

J o 

may be obtained from the formula given by Ng and Geller [4], p. 155, 3.5 (33), which may be put in the form 



F s (2P,a) = -^exp(-j3 2 ){erfi[j8] - erfi[/ta(a 2 + l)~^}. (3) 

For a < 0, one may use the relation 

erfc(-s) = 2 erfc(0) - erfcfc) (4) 

to obtain the Fourier sine transform connection 

F,(2j3, - \a |) = 2F s (2/3, 0) - F,(2j8, \a\). (5) 

The analogous Fourier consine transform, for a > 0, is defined here as 

F c (2/3, a) = ( cos(2i8.r)exp(-.r 2 )erfc(a.r) dx. (6) 

J o 

Two closed form special cases of this are 

F r (2l3.0) = (Vn/2)exp(-l3 2 ) (7) 

and 

F c (0,a) = (7r) _1/2 aictan(a- 1 ). (8) 

This latter equation follows from reference [1], p. 7, 4.3(2). We shall demonstrate that for the general case, 
the Fourier cosine transform may be written as 



an+a 2 r m 



(7r) _,/2 arctan(o:~ 1 ) + a I exp(a 2 x 2 ) erfi(x) dx 

J 



(9) 



F c (2(3,a) =exp(-)3 2 

or alternatively 

F c (2p,a) =a(7T)- m exp[-p 2 (l + a 2 )" 1 ] times f [z 2 + a 2 )~ l exp[i8 2 z 2 (l 4- a 2 )" 1 ] dz. (10) 

Jo 

These forms may be convenient for numerical evaluation when closed form expressions are not feasible. For 
Ot = 1, the integral in eq (9) may be evaluated (see reference [1], p. 7, 4.3(1)) to give the additional closed 
form special case 

F c (2p. 1) = (VW4)exp(-/3 2 ){l + 2[erfi(/3/V2)] 2 }- (U) 
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The general expressions given in eqs (9) and (10) may be established starting from the formula given by 
Erdelyi [2], p. 307: 



J x 2n cos 
J o 



(2 (3.x) erfc(ax) dx 



'> + ') .rY.4-1 -4.H..-4.S/. i/.. _£ 



2 F 2 | n + I, n + Va; n + 3 /2. Va; - - 
\/ir(2n+ ])a 2 " +l \ a 

1 ^ T(k + n + 1) „.,.*, 

= > (-6 2 /a 2 ) k 

2a 2n+1 n (ft + n + Va) r(ft + Va)*! P ' 

It then follows using the series expansion for exp(— x 2 ) that 

F (2/8 «) - 2 -t^i 2 - F(A+ ^ +1) (-/BVoV 

,IP ' ' n./!2a 2j+1 ft (A' + ./ + Vi)r(ft + V»)JB ' P 7 ' • 

Reversing the order of summation and employing the series representation of the hypergeometrie function, it 
follows that 

1 ^ (— B 2 /a 2 ) k 
F c (2j3, a) = — 2 n ; ' 3 A ^ + 1/2 - * + 1; * + 3 ^ " I A 2 )- 

ZQ! a-=o MA' T (2) 
Now, substituting the integral expression of reference [5] p. 114, 

| YL _|_ 3/ 2 ) r I 

2 F,(A- + l /2, A- + I; /: + 3 / 2 ; -l/a 2 ) - ' 1/x t*< I - /)" ,/2 ( I + t/a 2 )- k ~ l/2 dt, 

k\\ (V2) J () 

and summing the geometric series, one obtains 

b\.(2fi.a) = — ^= f ,^ 2 '<« 2 +'>-'(i - t )- ll2 (] + t?/a 2 )- ll2 dt. 

2a\/TT Jo 



From this, the formula quoted in eq (10) follows immediately upon change of the integration variable from / to 
z (using the positive square root branch forz) via 

(1 + t 2 /a 2 ) = (1 + a 2 )/(z 2 + a 2 ). 

The integral in eq (10), of the form 

r a (f) s f ^ 2(14fl2r V + a 2 )" 1 dz 
Jo 

satisfies the differential equation 



dg ( 1 + a 2 ) ,/2 

Hence, integrating this equation with respect to £ over the range (0, j8), recognizing that 7^(0) = a~ l 
arctan(a _1 ), and setting f = (1 + a 2 ) 1/2 x, one finds that the expression given in eq (10) may be written in the 
equivalent form displayed in eq (9). 

For negative a, analogously to eq (5), the Fourier cosine transform connection is 

F c (20, ~\a\) = 2F c {2p, 0) - F c (20, |o|). (12) 
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3. Evaluation of \(p, [3, 0; e, </>) 

The five parameter definite double integral defined in eq (1) will now be evaluated in terms of the Fourier 
transforms discussed in the previous section. Results are quoted in eqs (15), (17), and (18). The results for 
the special case )3 = are given in eq (19), while the results for the limiting case ^8 — > y as j8 — » are 
presented in eqs (21-22). 

We begin the derivation of the results by transforming the domain of the double integral using the general 
relation 

r. cc ~ oo ~ oo ~ y 

dj; I dr\ h{^ J]) = I dy I dx h{y — x, x), 

J o J o J() J 

which may easily be established using the variable changes f + r)^> y, 77— > :*;. Consequently, it follows that 

-00 r.y r.CC ~oo , 

I dy dx h(y\ x) = I dy I dx h(y + x, x). (13) 

J o J o Jo Jo 

Employing this relation, we may write /(p, )8, 0; e, (/)) in the form 

-.00 ,, 00 

2/(p, (3, 0; 6, </>) = J>m{\ dy dx e -^^>V^V^ + ^^^ ] 
J o Jo 

/. oo - oo 

Jo Jo 

At this point, /3 will be restricted to values greater than zero. [The special case /3 = will be discussed later.] 

In analogy with Jones and Klein [6, p. 3], we use a Dirac 8-function to write the double integrals as triple 
integrals. 

2/(p, /3, 0; e, <l>) = 

f 00 ,.00 /* 00 

^m{e iw+ * ) <fc e«" d r I d* ( -- p2( » + *>V* 2 8[t - >• - (1 + e)x] 

J - oo Jo Jo 

/» 00 ,. 00 ..00 

+ e*»-*> dt e m dy dx e~ pH y+x)2 e~ *g[t ~ y ~ (1 - e)x]}. 

J —oo J Jo 

Doing the integration over y first and noting that y > in the range of integration, we obtain 
2I(p, /3, 0; e, c/>) = 

3>m{e m ™ \ dt e m f dx e^^V^I) - (1 + e)x] 

J —cc Jo 

*. 00 ,.00 

+ e i,e -* ) dl e ,w dx e- p2< ' + " , V-- r2 0[< - (1 - e)x]} 

J —cc Jo 

where the step function 0(z) is defined as unity forz > and zero for z < 0. The presence of the step function 
restricts the relevant range of integration for the variables t and x in the following fashion. In the first of the 
preceding integrals, the step function requires that (1 + e)x < £, so since x must be > 0, it follows that t > 
0. Consequently, in that integral the ranges are 

< t < oo and < x < (1 + e) -1 *. 

In the second of the preceding integrals, the step function requires that (1 — e) x <t. Thus: 

(i) if e < 1, then < x < (1 - e)~H and < t < oo, 
(ii) if 6=1, then < t < oc and < a; < oo, 
(iii) if 6 > 1, then oc > x > max{-(e - \)~H, 0} and -co < % < oc. 

322 



It is useful to define the positive quantity 

IX =(] +eV)" 2 . (14) 

For convenience, in the first of the preceding integrals we change variables from (/, v) to (z, defined by 

z = (p/fJ,)t, f = fix - epz 

so that the appropriate range in the first integral becomes 

- ¥ 2 



< z < oo -epz > £ < 

Ld + e)pJ 

Similarly, in the second of the preceding integrals, variables are changed from (/. v) to (z, ^) now defined by 

z = (p/n)t, f = jUL.x + €/>Z 

so that appropriate ranges for the second integral for the various cases becomes 

1 + eg 2 " 



(i) if e < 1, then < z < », €/ > z < £ < 



( 1 - c)j> 



2, 



(ii) if e = 1, then < z < 00, epz < £ < 0°, 



(iii) if 6 > 1, then — oc < z < oc. 



+ ep 2 



(e - D/> 



z, e/;zf < £ < ». 



Making these substitutions in the preceding integrals we find 
2 P I(p,P,d; €,</)) = 



.r±v 



^00 - B(*) 

dz e im ' p)z e- z2 d£e-* 2 } 

JE J A(z) 



yher 



(i) if e < 1, then E = 0, A(z) = epz, £(z) = 



+ ep 2 



(1 - e)p_ 



z, 



(ii) if € = 1, then E = 0, A(z) = epz, B(z) = », 
(iii) if e > 1, then E = — », A{z) = max] — z, epzj^ , #(z) = w. 



(e- Dp 

From the definition of the erfc function and eq (4), it follows that 
(i) forO <€< 1, fi = (1 + eV) 1/2 : 



47r- ll2 pI(p,P,6; e, 0) 



^ 



+ 



2 { e i<*-w>) f 

X 

Jo 



dz e' w/p) V^ 2 



v 1 - ep 2 

2erfc(0) — erfc (epz) — erfc I — z 

(1 + e)p 



jiPMPte-z 1 



erfc (epz) — erfc I — z 

1 (1 - e)p 



h 
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(ii) for e= 1, fl= (1 + p 2 ) m : 

4ir- 1/2 pl(p,f3,0; 1, (/>) 



= ^TO 



^ gWM/P^ 



r 00 

+ e w ^> [ dz ^/^-* 2 erfc (pz)}; 
(iii) for e>l, /* = (1 + e 2 p 2 ) 1/2 : 



2erfc(0) - erfc(pz) - erfc ( — ^~ 



2p 



^m{e ue+<t>) ( 
Jo 



2erfc(0) — erfc(epz) — erfcl — z 

1 (1 +e)p 



| ( ' * 6p2) X 



»| dze i0lilp)z e-^er(c(€pz)}. 



Using the definitions of the Fourier transforms F c and F s given in the previous section, these results may be 
put in the form: 

(i) for <e< 1, fi = (l + e 2 p 2 ) m ; 
4<rr- ll2 pl{p,l3,0; e, <p) 



= sin(0 + </>) 
+ cos(6» + <f>) 



2F,.(pti/ P . 0) - Flpfi/p. 
2FJfin/p, 0) - FApii/p, 



1 - ep i 
(1 + «)fl 
1 - ep 2 
+ «)«/_ 



(15) 



- sin(0 - 4>)F e \BlL/p, T^rP) - cos{& - Wifin/P. £ + -Q 
+ 2sin(<p)[- cos(d)F c (Pn/p, ep) + sin(e)F s (fiix/p. ep)]. 

[Note: The Fourier transforms have been defined for non-negative values of their second arguments. When 
(1 — ep 2 ) < 0, the above expressions may be used by replacing 



'I Pn/p. ' ; €/ f ) = Flfiii/p, 1 



(1 +e)p/ V (1 +e)p, 

where F represents either the cosine or the sine transform. This replacement follows from eqs (5) and (12).] 
(ii) for e=l, fl = (1 +p 2 )" 2 ; 

47r 1/2 p/(p,/8, 0; i. <M 

i-/ 



(16) 



= sin(0 + <p) 2F c (PfJL/p, 0) - FJ /3/n/p, 

+ cos(0 + J 2F t (fipi/p, 0) - Flfiii/p, ^ L 

+ 2sin(<p)[-cos(0)F,O3p/p, p) + sm(0)F s (p[i/p, p)]. 
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(17) 



Again, for 1 — p < 0, the above expressions may be used if we make the replacement given in eq (16) with 
e= 1. 

(iii) for e> 1, /x = (1 + eV) 1/2 ; 

477- ,/ y (/'•£- 0; C </>) 

- e/> s 



= sin(0 + (/>) 2F c (Pfx/p, 0) - FJ £ M /p 



)Uf s 



4- cos(0 + (/>) 2F s (£ix/p, 0) - FJ j8/Lt//), 



(1 + e)fl 



"€/> 2 



(1 +€)/), 



18) 



/ 1 + ep 2 \ / 1 + e/r 

+ sin(0 - 0)fJ 0/i/p, — g-J - cos(0 - 0)F f l jB/i/p, ^7^- 

+ 2sin(<£) [- cos(0)F c (Pfx/p, ep) + sm(6)F 8 (pix/p, ep)]. 

Here, when 1 — ep 2 < 0, the replacement given in eq (16) should be made. 

The various results given for/(/>, /3, #; e (/)) are continuous for e — 1. Using the fact that F ( .{/3, x) — > and 
F s .(/3, x) — * as a — > x , we may readily verify that 

lira /(/>./3. 9; e. c/>) - /(/>.£. 0; I. 0) 

e— »1 

and 

lim /(p,/8, 0; 6, (/>) = l{/Kp. 0; L, (/)). 
For the special case (3 = 0, one obtains directly from the definition of the double integral under discussion, 

/(/>, 0, 0; €.(£>) = cfy ^vexp(-/> 2 v 2 )exp(-.v 2 )sinr^sinc/) 

r 00 

= sin#cos(/> r/yexp(— /; 2 y 2 )[eric(0) — eHc( v)] 

= Vtt (2p) _1 sin0cos<£ [F c (0, 0) - F c (0, I //>)]. 

Thus, it follows using the specific values of the transforms, 

I(p, 0, 6; e, 0) = (4p) _1 sin0 cos [77 - 2arctan(p)] 
= (2p)~ 1 sinO cose/) arctan (l/p). 



(19) 



This result also follows by a limiting process (3 — ■> in the results quoted in eqs (15-18) upon application of 
the relation for positive arguments 



and the identities (for e > 0, /> > 0): 

arctan (ep) ± arctan (/;) = ± arctan 
and 



arctan(x) + arctan(l/x) = 77/ 2, 
( I ± e)p 



arctan I — ) ± arctan 
P 



(iH 



arctan 



1 +/ 



(1 ± e)p 



ep< 



for 1 - ep 2 > 0, 



for ep 2 - 1 > 0. 
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[The inequivalent forms corresponding to the two upper signs are necessitated by the arctangent's definition 
which requires 77 /2 < arctan(x) < tt/2.] 
Finally, the particular double integral 

7/(y, p, <£) s I dy d*exp(-p 2 y 2 )exp(-* 2 )cos(yx + 0), y > (20) 

J >/ 

may most easily be evaluated by integration by parts: 

//(y, /;, 0) = —(1/p) I [derfc(y)] I d* exp(— * 2 ) cos(y:\: + </)) 

= Vk (2p)->[(co S <t>)F c (y. p) - (sincf>)F s (y, />)]. 
This result also follows directly from Eq (18) by the following limit process: 

H(y. p. </>) = lim Up. (3. tt/2; y/fi, </>). (22) 



(21) 
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